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ABSTRACT 
> ■ 

^ I Propagation of high-frequency radiation in an anisotropic medium can be 

I described in terms of two normal modes with different polarizations and different 

if^ \ absorption coefficients. We investigate the properties of the normal modes in 

a strongly magnetized hydrogen gas for conditions expected in atmospheres 
of isolated neutron stars. We use the Kramers-Kronig relations to obtain the 
polarizability tensor for the strongly magnetized hydrogen. We derive and 



O 

H i compute the polarizations and absorption coefficients of the normal modes from 



the polarizability tensor using both analytical approximations and numerical 
calculations. We find that the spectral features and anisotropy associated with 
the bound-bound and bound-free transitions in the magnetized hydrogen are 
^ [ manifested in the polarization characteristics, which affects substantially the 

spectral and angular dependences of the absorption coefficients of the normal 
modes and the transfer of radiation in neutron star atmospheres. There exist at 
least two critical frequencies where, for any direction of propagation, either the 
normal mode polarizations are exactly linear or orientations of two polarization 
ellipses coincide with each other. For each of these frequencies, there exist a 
direction of propagation for which the two normal modes are linearly polarized, 
and their polarizations and absorption coefficients fully coincide. The unusual 
properties of the normal modes should manifest themselves in the spectra, 
angular distribution, and polarization of the thermal-like radiation emitted from 
surface layers of neutron stars. 



^On leave from A. F. loffe Physico- Technical Institute, St. Petersburg, 194021, Russia 
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1. INTRODUCTION 

With the development of observational capabilities of space X-ray/UV observatories 
and ground based telescopes it is now feasible to detect radiation from surface layers of 
neutron stars with rather low temperatures, T ~ (0.2 — 10) x 10^ K (Ogelman 1995; 
Pavlov, Stringfcllow & Cordova 1995). The immense magnetic fields of neutron stars, 
B ~ 10^^ — 10^'^ G, substantially increase the ionization potentials of atoms (e. g., Canuto 
& Ventura 1977): / ~ Z^Ry In^ {-f / Z^) at 7 > Z^, where 7 = huBe/ i2Ry)= 5/(2.35 x 10^ 
G) is the magnetic field in the atomic units, = eB/rrieC is the electron cyclotron 
frequency, Z is the effective ion charge, and Ry = mee^/2/i^ = 13.6 eV is the Rydberg 
energy. This means that even light atoms may be only partly ionized (or even fully 
nonionized) in neutron star atmospheres, and be the main contributors to the atmosphere 
opacity (Pavlov et al. 1995a). Thus, to model the atmospheres for further interpretation of 
the observational results, it is necessary to investigate the effect of the strongly magnetized 
atoms on absorption and propagation of radiation. 

It is well known (see, e.g., Gnedin & Pavlov 1974) that radiation propagates in an 
anisotropic medium in form of two so-called normal modes (NMs) which have different 
polarizations and, consequently, different absorption and refraction coefficients. The 
dependence of the NM absorption coefficients on the frequency and direction of propagation 
determines the spectrum and angular distribution of radiation emitted from an anisotropic 
medium. Properties of the NMs have been thoroughly investigated (e.g., Pavlov, Shibanov 
& Yakovlev 1980; Kaminker, Pavlov & Shibanov 1982, 1983; Meszaros 1992, and references 
therein) for the case of fully ionized plasma, when polarization and absorption of the 
modes are determined by the free-free transitions and Thomson scattering in the magnetic 
field. However, there has been no analysis of how these properties would change if the 
bound-bound and bound- free transitions in the strongly magnetized atoms are important. 

To find the polarizations of the normal modes, one should proceed from the 
polarizability tensor of the anisotropic medium. The contribution of the bound-bound and 
bound-free transitions to the components of the tensor is determined by the structure of 
the strongly magnetized atoms. For the simplest hydrogen atom, assuming the center of 
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mass of the atom is fixed, the atomic structure and probabihties of the radiative transitions, 
which determine the antihermitian part of the polarizabihty tensor, have been investigated 
in detail (e. g., Canuto & Ventura 1977, Forster et al. 1984, Potekhin & Pavlov 1993). The 
hermitian part of the polarizabihty has not been discussed, to the best of our knowledge, 
even under this oversimplified assumption. Moreover, it has been understood recently that 
motion of atoms in strong magnetic fields changes qualitatively their structure (Herold et 
al. 1981; Vincke & Baye 1988; Vincke, Le Dourneuf & Baye 1992; Pavlov & Meszaros 1993; 
Potekhin 1994). This occurs because the motion across the magnetic field induces an electric 
field, F = c~^V X B, which breaks the cyhndrical symmetry and couples the center-of-mass 
motion to the internal structure of the atom. The energies and wave functions of the moving 
atom depend on the transverse component of a generalized momentum K. When 
is small (at low temperatures of the gas), the perturbation approach is applicable, and 
the coupling effect can be described in terms of an effective transverse mass of the atom 
(Vincke & Baye 1988; Pavlov & Meszaros 1993). This transverse mass always exceeds 
the mass of the atom without the magnetic field. It grows with the magnetic field, being 
and is generally higher for more excited atomic states. When increases beyond the 
applicability of the perturbation theory, the induced electric field causes decentering of the 
atom, i. e., the electron gets into a potential well which is shifted from the proton more for 
the larger K^. These dcccntcred atoms have very unusual properties. For instance, their 
velocity decreases with increasing Kj_ because the growth of K_|_ goes not to acceleration 
of the atom but to larger decentering. The decentered states have relatively small binding 
energies, so that they can be occupied only at relatively high temperatures and not too 
high densities, when the size of the atom does not exceed the mean distance between the 
particles. 

Qualitative effects of the coupling on the radiative transitions were discussed by Pavlov 
& Meszaros (1993). They show, in particular, that different dependence of different atomic 
levels on leads to a new broadening mechanism of the spectral lines and photoionization 
edges. The "magnetic coupling width" exceeds the Doppler width by orders of magnitude, 
growing proportional to the gas temperature. More detailed treatment of the radiative 
transitions, with allowance for the decentering of moving atoms, has been presented by 
Bezchastnov & Potekhin (1994) for bound- free transitions and by Pavlov & Potekhin (1995) 
for bound-bound transitions. 

In the present paper we will consider the polarizabihty tensor and properties of the 
polarization modes of the hydrogen gas for the case of relatively low temperatures, when 
interaction of the radiation with the medium is mainly determined by the nonionized atoms. 
At these low temperatures, most of the atoms occupy the ground level, and, at densities 
typical for neutron star atmospheres, one can neglect contribution of the decentered states. 
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In these conditions, the couphng effects can be considered in the frame of the perturbation 
approach. In Section 2 we recall general equations for the polarizability tensor, normal 
mode polarizations and absorption coefficients, and describe our approach for calculation 
of these quantities. We present numerical results and analytical limits in Section 3. In 
Section 4 we discuss some consequences of the obtained results for neutron star atmosphere 
modeling. 



2. BASIC EQUATIONS 



2.1. Polarization Characteristics of the Normal Modes. 



The propagation of electromagnetic waves in a magnetized medium is fully described 
by the dielectric tensor 

Qfc = 8ik + 47rxifc , (1) 

where Xik is the polarizability tensor. In the case of a tenuous medium, (i. e., at \nj — 1| < 1, 
where j = 1,2 labels the NM, and Uj is the complex refraction index) the NMs are 
transverse, and the (real) refraction index Kj and the absorption coefficient are given by 
(Pavlov et al. 1980) 

= + i—iij ^ni± ^/nl + ri^ , (2) 

where 



71/ = 1 + 71 {xyy + Xxx cos^ 6 - Xxz sm 29 + Xzz sin^ 6*) , (3) 



riL = T^iXyy - Xxx cos^ 9 + Xxz sin 29 - Xzz sin^ 9) , (4) 

nc = 2ni {xxy cos 9 + XyzSin9) . (5) 

Here 9 is the angle between the wave vector q and the magnetic field B, and Xik is the 
polarizability tensor in the coordinate frame with the z-axis along B and the x-axis in the 
B-q plane. The NMs can be described by the ellipticity Vj (the modulus \Vj\ is the ratio 
of the minor axis to the major axis of the polarization ellipse; the sign of Vj determines 
the direction of rotation of the polarization vector) and the position angle 6j between the 
major axis of the polarization ellipse and the projection of B onto the plane perpendicular 
to q. The ellipticity and position angle are determined by the modulus Vj and phase ipj of 
the complex quantity 

-nc±^nl + nl 
Tj eK^iiLpj) = ; (6) 
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r, + l' ^ 2 ■ 

It follows from these equations that the NM properties can be conveniently analyzed with 
the use of real parameters q and p defined as (Gnedin & Pavlov 1974) 

q + ip ^ — . (8) 
nc 

The NM polarizations are nearly orthogonal ~ 5i ± 7r/2; 5j ~ or 7r/2; and V2 — —Vi) 
when 2\p\ <^ + + 1. In particular, they are orthogonal when p^ <^ q^ (the most common 
case, which is usually realized in a wide range of frequencies and angles), or when p^ ^ q^ 
and \p\ ^ 1. In the case when p^ <C g^, the (orthogonal) NMs are linearly polarized, 
YPi\ = 1^2! — 0, if \q\ ^ 1, and they are circularly polarized, {Vil = {7^21 — 1, if \q\ <^ 1- If 
S> q^ and q^ -C 1, then the NMs are linearly polarized at \p\ > 1, and the polarization 
ellipses coincide [Si ~ (^2 ~ (7r/4) sign(p), Vi = —V2] at \p\ < 1. In the special case when 
g = and \p\ = 1, the NM polarizations are purely linear and coincide with each other (are 
"completely non-orthogonal"): Vi — V2 — 0, 61 — 82 — (7r/4) sign(p). For further details 
see, e. g., Pavlov et al. (1980). 



2.2. General Equations for the Polarizability Tensor in the Dipole 

Approximation. 



In the dipole approximation, which is justified for frequencies that are not too high (see, 
c. g., Pavlov & Potckhin 1995 for a discussion), the polarizability tensor of a magnetized 
gas is diagonal in the cyclic coordinates, Xai3 = XaSa/3, where a,P = 0, ±1, and the cyclic 
unit vectors are defined as eo = e^, e±i = 2~^/^(ej; ±iey). In terms of the cyclic components 
Xa, we have 

l + cos^6l , . 2 



rij — 1 + TT 



ix+i + X-i) 



+ Xosin 9 



( X+i + X-i 
V 2 



- Xo sin^ 9 , 



and 



q + ip^ {q + ip) 



n-c = 7r(x-i - x+i)cos6l , 
sin^^ 



2cos^ 



X+i + X- 
X-1 - X+1 



2X0 



(9) 

(10) 
(11) 

(12) 
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The polarizabihty tensor is the sum of the hermitian and antihermitian parts, Xa = Xa > 
responsible for refraction and absorption of radiation. The antihermitian part can be 
written (for a single component ideal g as) as 



Xa 



UJ. 



^Ee. frfS{uj-L0f,)-ff,6{uj-uj,f) 
?^E'(e. - 0) \JTf - ^n) - f% - ^^f) 



UJ, 



6uj 



(13) 
(14) 



if 



where ujpa = (47rnae^/me)^/^ is the "plasma frequency" of bound electrons, is the number 
density of atoms, = Uai/ua is the population of the quantum state i, ujji = {Ej — Ei)/Ti is 
the transition frequency, and fij is the oscillator strength for the transition i —>■ f, 



fi 



2meUJfi 



-fi 



■ en 



2meUjfi 

n 



|(x. 



■Oi)fi\ 



(15) 



The double sum J2if implies summation (and integration over the continuum part of the 
energy spectrum) over all possible combinations of i and / (except i = f); in the sum X^i/ 
only the states are summed (integrated) for which Ef > Ei. The general equation for the 
hermitian part can be written, in the same approximation, as 



X, 



H 



UJ. 



pa 



if 



fa 
Jif 



f 



fi 



UJ — UJfi UJ — UJif 



UJ. 



pa 



iTlUJ 



E'(e. - 0) 



if 



fa 
Jif 



f 



fi 



UJ - UJfi 



UJ 



UJ, 



if 



(16) 



It follows from the completeness of the wave functions that 



E- /5 + //" 

/ ^ft 



0, 



(17) 



which means, in particular, that Xa has no pole at = if all the electrons are bound. 

We see from these equations that direct calculation of x^ implies an additional 
integration over the continuum states, as compared to x^? because there are no delta- 
functions in equation (0), and this integration may not be straightforward. An alternative 
way to find the hermitian part is to use the Kramers-Kronig (KK) relation (Guilbert 
transformation of the antihermitian part): 



Xaiuj) = - 
TT 



xii^') 



oc, UJ — UJ 



duj' ^ Kixiiuj)] 



which is derived from the general analytical properties of the polarizabihty: x{^) is a 
regular function in the upper half-plane of complex oj; x{^) has no singularities on the 
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real axis except possibly a simple pole at the origin; x(c<j oo) — > (see e.g. Landau & 
Lifsliitz, 1960). Since the quantity uJXa{^) also has all the necessary analytical properties, 
we may apply the KK transformation to this product 

(co) = -Ki^xti^)] = T^Klf-^aico) e(c.) + 0(-c.)] , (19) 

where 9(x) = 1 for x > 0, 6(x) = for x < 0, Haiuj) is the absorption coefficient (defined 
for u > 0): 

2 

/i„(u;) = - 0) /f/ '^(^ - ^/O = ^41 - exp(-^a;/A;T)] 5:e.o^.(^) , (20) 

if i 

where ai{uj) is the absorption cross section, and T is the temperature. To derive equation 
(p!9|), we used the property of the oscillator strength 

ftf = -fjr (21) 

which implies the symmetry relations 

Xt{-^) = -XU^), (-^) = X^J^) • (22) 

Applying the KK transformation to cJXa instead of Xa is technically easier because the 
additional u cancels the false pole (for bound electrons) of Xq, at w = 0. Thus, to find x^ in 
this way, one needs the absorption coefficients at u; > for the longitudinal (a = 0), right 
(a = +1), and left (a = — 1) polarizations. 



2.3. Oscillator Strengths and Cross Sections for a Strongly Magnetized 

Hydrogen Atom at Rest 

The absorption coefficients for the bound atoms are contributed from the bound-bound 
and bound-free transitions. For isolated atoms at rest, the problem of calculating the 
bound-bound absorption coefficients is reduced to calculating the oscillator strengths /"j. 
When the magnetic field B = Be^ is so strong that the ratio of the electron cyclotron energy, 
hujBe = 27Ry = 11.6i?i2 keV, substantially exceeds the binding energy, ~ Ry ln^7 at 
ln^7 ^ 1, the adiabatic approximation is valid. In this approximation, the bound states of 
the hydrogen atom are described by the Landau number = 0, 1, 2, the negative of the 
2;-projection of the angular momentum s = —N, —N + 1, — A^ + 2, and the longitudinal 
quantum number n = 0, 1, 2, ... (see, e. g., Canuto and Ventura 1977 for details). The wave 
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functions are the products of the Landau states $Ars(r±) and longitudinal wave functions 
QNsniz) which obey a one-dimensional Schrodinger equation. The oscillator strengths of 
the transitions between the states k — {Nsn) and k' — {N's'n') are given by the following 
equations 

= ^[(iV + .)W+(iV+l)5^,^+i] (23) 
Z^' = ^[iN + s + l)5^,^ + N5^,,r^^,]\z;^l\%>,s+,, (24) 
fL' - ^S,J^'N\zil\\ (25) 



where 



dz g*Az)g,{z) and = / dz g*Az){z/aB)g.iz) (26) 

-OO J —OO 



are the overlapping integral and dimensionless matrix element of the longitudinal coordinate, 
respectively; as is the Bohr radius. The transition frequency is 

oJk'k = (N' - N)ujBe + (N' -N + s'- s)u;bp + Cjk'k , (27) 

where huj^'n = ^k' ~ is the difference of the longitudinal energies, and ujbp = {iTie/i^p)^Be 
is the proton cyclotron frequency. In the application to the neutron star atmospheres, the 
only important case is when the initial state k belongs to the ground Landau level N — 
(the excited Landau levels are unpopulated at kT <^ hcuBe)- Moreover, we will consider 
temperatures low enough that the population of even the first excited level (010) is small 
compared to the ground state, (000), which allows us to neglect transitions between the 
excited levels. Then, putting k = (000), we obtain 

f-l CUge + ^l-ln',000 , ^(0) |2r r /r,oX 

Jkk' - ^ \ l^l-ln',000l CiN'\Os'-\ , [^o) 



Jkk' — ^ l^01n',000l ^iV'O^^s'l , \^^) 



rO _ ^^OOn',000 I ^(1) |2r c /onN 
~ — Ry — I 00n',000l On'oOs'O ■ WUJ 

Note that the longitudinal energies and wave functions of the states (Oln) and (1-ln) 
coincide (Hasegawa & Howard 1961). This means that a;i_in',ooo = '^oin'.ooo and 
■^i-in'.ooo = -^oii',000- Since the parity of the longitudinal states is (—1)", transitions 
from the ground (even) state are allowed to odd states, n — 1,3,..., for the longitudinal 
polarization a — 0, and to even states, n = 0, 2, . . ., for the circular polarizations a — ±1. 
Moreover, in the strong magnetic field, the overlapping integrals -^oi^j/^ooo n' are 
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much smaller than that for n' = and can be neglected. As a result, equations (p8[) - ( pOD 
take the form 

fi^K' = 2(1 + j 6N'lSs'-lSn'0 , (31) 

.+1 ^ / '^010,000 , ""^e] . . . ,„r,\ 



UJBe 



V, 



~ — Ry \^oon',ooo\ ^N'o(>s'o ... {66) 

These equations demonstrate that the bound-bound transitions are quite different for 
different polarizations. The strongest transition (with ~ 2, oj^'k, — ujse) occurs onto 
the excited Landau level, (000) (1-10), and it is allowed only for the left polarization. If 
the electrons were not bound to protons by the Coulomb forces, — >■ 0, this transition 
would turn into the usual electron cyclotron absorption. Although the absorption in the 
left polarization is negligible for the most interesting range uj <^ ujbb, we shall see below 
that it plays a crucial role for correct evaluation of x±i and of the NM polarizations in the 
low-frequency range. The only important transition for the right polarization occurs onto the 
first excited tightly bound level (010); its oscillator strength, ~ '2['y'^ ln7 + (mg/mp)], 
and transition energy, ^ci^^'k ~ Ryln7 + huBp, are much smaller than those for a = —1. 
This transition would be the proton cyclotron absorption if Coulomb binding could be 
neglected. The transitions allowed for the longitudinal polarization proceed onto the 
hydrogen-like levels of odd parity, creating a series of spectral lines with the leading line 
formed by the transition to the lowest hydrogen-like level, (000) (001). The oscillator 
strengths, /^k' < 1, and the transition frequencies, /icu^'K ~ ^Y^t^I, are smaller than for 
a = —1 but substantially greater than for a = +1. These transitions, solely caused by the 
Coulomb binding, disappear when the binding is turned off (if the electron-ion collisions are 
neglected) . 



Equations (PTD-(|55D are obtained in the adiabatic approximation, which leads to some 
errors in the oscillator strengths. In particular, this approximation implies that cj^'k ^ ^^Se 



so that one cannot trust second term in the numerator of (|3T|) . Indeed, the slightly different 
approach (still in the frames of the adiabatic approximation) used by Hasegawa & Howard 
(1961) leads to the minus sign before this term, instead of plus. Moreover, equations 
(0) and (0) do not exactly satisfy the sum rule for the circular polarizations (see, e. g., 
Hasegawa & Howard 1961), 



E [f^.' + = 2 . (34) 

k' 

Indeed, the sum of the oscillator strengths (PTT) and (|5^ for the two transitions, even 



without account for, e. g., bound- free transitions for the circular polarizations (see 
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below), exceeds 2. Therefore, more accurate values for the oscillator strengths, obtained 
numerically with allowance for the non-adiabatic corrections, should be used for precise 
calculations. We shall use below the values from Forster et al. (1984). For the magnetic 
field B = 2.35 x 10^^ G (7 = 1000) and for the transitions included into our calculations, 
they are listed in Table 1. 

For the bound-free transitions, the final state belongs to the continuum: k' = {N's'e'u'), 
where e' (> 0) is the electron energy in the (one-dimensional) continuum, and i^' = ± is the 
parity of the longitudinal function. The sum over the final states in equation (EO) is 



^ = ^2^ 5: 77^7' (35) 



N's'u' 

where is the the z-extension of the periodicity length (the normalization length) of the 
continuum wave function of the final state. The photoionization cross section is the sum of 
the partial cross sections 

a:{uj) = <^N's'A^) (36) 

N's'u' 

which differ from zero for frequencies above the photoionization thresholds, 
o;^]^, = eooo/^ + [^'(1 + ^e/^p) + s'me/mj\ujBe if only the ground state is populated. The 
photoionization cross sections for the atoms at rest have been considered by Potekhin and 
Pavlov (1993) whose approach and codes we shall use for computations and analysis. In the 
dipole approximation, only one term contributes to the sum (|36|) for a given polarization: 
N' = 0, s' = 0, z/' = - for a = 0; A^' = 1, s' = -1, z/' = + for a = -1; and A^' = 0, s' = 1, 
z/' = + for a = +1. Then the cross sections can be expressed in a form similar to that for 
the bound-bound transitions: 





2 

= Tia^ap 


1 1 
l^l-le'+,000l 




2 ^CU 

= na^ap 


17(0) |2 
I^Ole'+.OOOl 


^» 


2 1 
= Tianap — 

^ ^Ry' 


7(1) 12 
^00e'-,000l * 



Ry 



Ry 



— , (38) 

(39) 



ON ' 



where ai? is the fine structure constant. 



^Vs'eV',000 = {Ljas) / dzglj,,,^,^,{z) {z/aefgomiz) {k = 0, 1) , (40) 

J —00 

htu^ = eooo + ^t^Se, ^^^t^ = eooo + ^t^Sp, and hcu^ = eooo are the threshold energies, and 
e' = h{uj — tuf ) in the continuum wave functions QN's't'u'iz)- Similar to the bound-bound 
transitions, we have ^i^^^v 000 ~ ^oievooo- shown by Hasegawa & Howard (1961) 
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and Schmitt et al. (1981) that the overlapping integrals Z^^^ and the matrix elements Z^^^ 
are proportional to (cu — uj^y^"^ at cu — >^ a;" + 0, so that the cross sections are finite near the 
thresholds. The overlapping integrals, and the cross sections for the circular polarizations, 
are very small because the continuum wave functions are almost orthogonal to the wave 
function of the ground state. A convenient measure for the strength of the photoionization 
absorption for a given polarization is the cumulative (integral) bound-free oscillator strength 
obtained by integration over the final continuum states. Prom Table 1, where the integral 
bound-free oscillator strengths are presented together with the photoionization thresholds, 
we see that the bound-free transitions are much weaker than the bound-bound ones for 
the circular polarizations, being somewhat stronger than the bound-bound transitions 
for the longitudinal polarization. Thus, the main transitions arc: (000) — > (OOe'— ) and 
(000) (001) for a = 0, (000) ^(1-10) for a = -1, and (000) (010) for a = +1. Wc 
include all the transitions listed in Tabic 1 into our calculations of the absorption coefficients 
and the components of the polarization tensor. Note that the oscillator strengths in Table 
1 satisfy the dipole sum rules with accuracy certainly better than 1%. 



2.4. Broadening Mechanisms and Effects of Finite Atomic Velocities 



The above equations correspond to the case of isolated atoms at rest. They yield 
infinitely sharp features in both the antihermitian and hermitian parts of the polarizability. 
Thermal motion of atoms and interaction with surrounding particles broaden and smooth 
the sharp features. It was shown by Pavlov and Mcszaros (1993) that the most important 
broadening mechanism is associated with different coupling of the internal and centcr-of- 
mass motions (the difference between transverse masses, in the perturbation approach) of 
different levels. In particular, an infinitely narrow bound-bound profile 5 (a; — uJk'k) turns, 
in the perturbation approach, into 

r/ \ I r \ i (rM)"^exp[a; -u;«/«)/rM] a; < , . 



The magnetic width Tm is 



M 




(42) 



where is the transverse mass of the level k. The magnetic width grows hnearly with 
temperature {TlTm — kT for the transitions from the ground level to highly excited states) 
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and dominates all other broadening mechanisms. The magnetic broadening also smoothes 
the photoionization edges: 



/x^'(^) = na[l - eM-h^lkT)] / du;>M(^ - , (43) 

Jmax{w,aj°) 

SO that the bound-free profiles acquire exponentially decreasing low-frequency tails. 



Although the magnetic width is very large, we see from equation (|4T| ) that it still 
retains infinitely sharp high-frequency edges of the bound-bound profiles. In order to avoid 
logarithmic singularities in x^, arising after the KK transformation of such functions, one 
should take into account additional broadening mechanisms: Doppler broadening and/or 
collisional broadening. 

The Doppler width of the spectral lines becomes anisotropic in the presence of strong 
magnetic field (Pavlov and Meszaros 1993). For a transition k k', the Doppler width is 



nie kT 



where M = rrip + irte is the usual "longitudinal" mass of the atom, and 9 is the angle 
between the magnetic field and the wave vector. 

Since the density in the neutron star atmospheres is rather high, collisional broadening 
may be the most significant. The theory of collisional broadening in very strong magnetic 
fields has not yet been developed. A particular case of the Stark broadening by plasma 
electrons, which is expected to be the most important in strongly ionized gases, has been 
considered by Pavlov & Potekhin (1995) in the impact approximation. In the case of 
weakly ionized or non-ionized gases the broadening is mainly caused by interactions of 
neutral particles (e. g., van der Waals interactions). We know from the theory of spectral 
lines in the non-magnetic case (e. g., Sobelman, Vainshtein & Yukov 1981) that the 
most important broadening mechanism for strong transitions in a single-component gas is 
so-called self-broadening (or resonant broadening). It results in a Lorentz profile with a 
width that can be estimated as 

Ts ~ Jkk' TLa = 4- f^^> [naaB) ^ , (45) 

where is the oscillator strength of the transition, and is the number density of 
atoms. We adopt the same estimate for strongly magnetized atoms. The self-broadening 
width becomes negligibly small for very weak transitions. However, it seems obvious that 
the collisional width cannot be smaller than = F^^ + F^^' , where 

4 n (av) Ry 

~ 2na{av} ~ 2na(o-_LK(w±K) + (y\\^,{v\\)) = — ("-aOij)^ ^ , (46) 



-13- 



and fTj^K ~ 2(rK + ro){lK + Iq) and (7||« ~ 7r(rK + Tq)^ are "geometrical" transverse and 
longitudinal cross sections for a collision of the atom in the quantum state k with a 
background atom in the ground state. The cylindrical radius of the atom can be estimated 
as ~ aB[(|s| + l)/7]^/^; its longitudinal extension is 1^ ~ 1.2aB(|s| + l)^/^(ln7)~^ for the 
tightly bound levels, or ~ 1.6aB{lnt[(n + l)/2]}^ for the hydrogen-like levels (Potekhin 
1994). The transverse and longitudinal mean relative velocities of the colliding atoms 
are {v^^) = (7rA;T/2/x«)V2 and {v\\) = 2{kT /TiMf/'^ , where /x« = MlMl/{Ml + Ml) is 
the reduced mass. In addition, the lines are broadened by radiative damping (natural 
broadening). The corresponding widths are negligible for the transitions within the Landau 
level N — 0, but the radiative width considerably exceeds the coUisional widths for the 
main transition in the left polarization, (000)^(0-10). For this transition we have 



8 3 2% 



(47) 



The coUisional and radiative broadening of an infinitely narrow line results in the Lorentz 
profile 



27r 



(a;-a;«;/«) + 



(48) 



An example of the values of the magnetic, Lorentz and Doppler widths is presented in 
Table 1 for the transitions of interest. We see that the Doppler and Lorentz widths are 
generally comparable. Nevertheless, we shall neglect Doppler broadening in our calculations. 
The convolution of the Doppler and Lorentz profiles coincides with the Lorentz profile far 
from the line center (at \u! — cu^'i^l > Tr,), and the resolution in the calculations presented 
here is not sufficient to describe the details near the fine cores. Moreover, the above 
estimates of the widths are very approximate, and the primary object of introducing the 
broadening is to avoid sharp edges in the absorption coefficients. Thus, we assume that the 
bound-bound lines are broadened by magnetic and Lorentz broadening: 



nieC 




exp —jj^ 



/ y J K 



da;''0M('^ - - (^k'k) ■ (49) 



3. RESULTS 



3.1. Components of the Polarizability Tensor 



As a representative example, we investigated the NM properties for conditions 
typical for photospheres of moderately old, isolated neutron stars (Pavlov et al. 1995ab): 
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B = 2.35 X 10^2 G (7 = 1000, huBe = 2000 Ry, huBp = 1.09 Ry), fcT = 1 Ry 
(T = 1.58 X 10^ K), p = 1 g cm~'^ {naCi% = 0.0885). The ionization energy of hydrogen 
in this field is 15.32 Ry. All the relevant bound-bound and bound-free transitions, the 
oscillator strengths, and the widths are listed in Table |I[ Although the problem of ionization 
equilibrium in strongly magnetized, dense hydrogen has not been satisfactorily solved so far, 
preliminary estimates show that the ionization degree is expected to be low at this relatively 
low temperature and high density, which allows us to neglect the ionized component and to 
demonstrate qualitative features of the NMs in a nonionized, strongly magnetized gas. 

Three absorption coefficients /ia(c<j), calculated as described in Sections 2 and 3, are 
shown in Figure 1, and the corresponding antihermitian components of polarizability tensor 
Xai^) = (c/Attu) fia{^) are shown in Figure 2. We see that the behavior of and Xai^) 
is quite different for different a. For a = — 1, it has one very strong resonance associated 
with the bound-bound transition (000)^(1-10) and centered very close to the electron 
cyclotron frequency. The strongest absorption in the most interesting frequency range 
cij ^ dj^e is that for radiation polarized along the field {a = 0). It appears as a broad 
resonance around ~ 15 Ry with a fine structure near the maximum associated with the 
peaks of the overlapped bound-bound and bound-free absorption. The high-energy tail of 
the resonance is due the photoionization (000) (OOe'— ), and the low-energy tail is caused 
by the magnetic broadening of the bound-bound transitions [mainly, of (000) (001)]. 
Two features are seen in /i+i(c<j) and x+ii^)'- the broadened, asymmetric spectral line 
(000) (010) and the (much weaker) photoionization spectrum (000) (01e'+). 

If the gas were fully ionized, the X-i component would be almost the same — the 
peak would shift from 2004 Ry to the electron cyclotron resonance at 2000 Ry. The x+i 
component would turn into a narrow peak at the proton cyclotron frequency, the width 
of the resonance being determined by the proton-proton collisions (free charged particles 
are not subject to magnetic broadening — their "transverse masses" are infinite). The 
longitudinal component Xo would decrease monotonically with u (the only "resonance" 
lies at = for the unbound particles), describing the inverse bremsstrahlung due to the 
electron-ion collisions. Thus, the Coulomb binding essentially alters the antihermitian part 
of the polarization tensor, especially for the longitudinal component. 

We made the KK transformation of c^Xa (see eq. [0) and obtained the hermitian 
components Xa (^) which are shown in Figure 2 for u > 0. The accuracy of the numerical 
KK transformation was tested with a simple model of fully ionized plasma. As expected, 
the hermitian part has spectral features associated with those of the antihermitian 
part. Indeed, as follows from equations ([ISD and ([19D , an absorption feature of Xai^) 
centered at a positive frequency uja, generates the contribution to Xai^) equal to 
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^paf°'[^'^^a{^a — ^)]~^ if the distance \u! — Ual from the feature is much larger than its 
effective width, and exp{—huJa/kT) <^ 1. Each spectral feature at the positive frequency 
uja has its "counterpart" at the negative frequency —u^a] e. g., for the spectral line 
(000)— ^>(010) in the right polarization {to+i = u;oio,ooo5 f~^^ = /ooo,oio) the counterpart 
line is centered at —u^i = — co'i_io,ooO) and its oscillator strength is = /ooo,i-io- This 
counterpart generates Ljp^f^"'[8'KUj_a{^^-a + ^)]~^- Summing these two terms with allowance 



for the property Uaf ° = oj-af^ (see eqs. [|T7l and l2TI]), we obtain the total contribution of 
the feature to the hermitian part: 

Svr {uj^-uj){uj_^ + u;) ^ ' 

This expression generalizes the well-known result (e. g., Davydov 1976) of the isotropic case, 
when oj-i = = Uq and = f^^ = /" are the frequency and the oscillator strength of 
a given transition. Notice that equation (|50D describes also the polarizability of a "cold" , 
magnetoactive plasma (e. g., Ginzburg 1970) if cUpa = cOpe = {Anriee^ /rrieY^'^ is the electron 
plasma frequency, tu^i = use, i^+i = ^^Bp, "^o = 0, /^^ = 2, f^^ = 2me/mp, and /° = 1. In 
our case of the cold, nonionized gas, x±i is determined by the bound-bound transitions, so 
that equation ( ^0]) immediately gives the circular components of the hermitian part: 

" 4.r(^_±"j)V,T^) ' ^''^ 

where uj+ = a;oio,ooo and U- = ci;i_io,ooo- The only difference from the corresponding 
components for the fully ionized plasma is that we replace by u^, and use by uj-. 

Interaction of the longitudinally polarized radiation with the cold gas is mainly 
determined by the spectral line (OOO)-^(OOl) and the photoionization (000)— > (OOe— ). The 
former yields 

H(bb) _ ^pafo 

where /o = /oooooi = ^^ooi.ooo- Accurate values for the bound- free contribution can 

be obtained only numerically. A qualitatively correct behavior of Xo'^'' can be obtained from 
an approximate model cross section: 

a\u) ^(Tt {ut/uf Q{u~ut) , (53) 

where Ut is the threshold frequency and ctj is the threshold cross section. Substituting this 



model expression into equations ( p!9D and (|20|) , we obtain 

^ , (54) 
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where f^^^^ is the integral oscillator strength of the bound-free absorption. At u <^ ut 
the polarizability resembles that for the bound-bound transition with cuq (3/2) ^/^cuj 

r - (3/4)/(°M) : 

3^^pa/(bf) 



Xo 



167r(|u;2 



When uj approaches Ut, the polarizability diverges logarithmically: 



^pa/(bf) 



In- 



2\iOt — Lo\ 



(55) 



(56) 



At C(j > cut, it first decreases, changes its sign at ~ 1.16co't, reaches a minimum at 
a; ~ 1.58a;t, and approaches zero at ljj: 



^pa/(bf) 



47r[u;|ln(cuVcuf) 



(57) 



Equations (|5^) through (^) are not applicable too close to ujj^ (=5.14 Ry), ujq 
(=14.4 Ry), ujt (=15.3 Ry), and uj- (=2004 Ry), where the bound-bound or bound-free 
polarizabilities diverge if broadening mechanisms are not included. The magnetically 
broadened spectral line generates the following longitudinal component of the hermitian 
part of the polarizability (at exp{—hujQ/kT) <^ 1) 
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(58) 



where Ei(a;) is the exponential-integral function. At \ujq 
equation (|52| ) with a shifted central frequency: 
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^paJO 
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^ Fm this equation turns into 



(59) 
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UJ^ 



Near the central frequency of the unbroadened line, the polarizability diverges 
logarithmically: 
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(60) 



StCUJqTm I \UJo—UJ\ 

where C ~ 0.577 is the Euler constant. Thus, with frequency growing from zero, the 
polarizability generated by a magnetically broadened, isolated spectral line first grows, 
reaches its maximum at ~ c^o — 1.35rjv/, changes its sign at u; ~ — 0.37rAf, reaches its 
logarithmically infinite minimum at uj = ujq, and then grows steadily to zero. 
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The magnetically broadened circular components are given by similar equations, 

,2 



pa 



iTCUJ 



exp 



Ei 



M± 



M± 



U0± — UO 
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Ei 



iOzfz + iO 



(61) 

which turn into equations ( pT]) far from the frequencies uj±. In the most interesting region 
u -C uJBe, the left-circular component varies very smoothly, 

(62) 



00. 



pa 



47rco'_(u;+ + oo) 

whereas the right-circular component shows a resonant behavior near oo^, similar to the 
above- described behavior of the longitudinal component around ooq. 

Magnetic broadening of the bound-free polarizability does not change the polarizability 
at \oo — oot\ ^ T M (see eqs. ^1[], and [0), and it removes the logarithmic divergence 
at oot: 



2 sO 



at 



oo 



OOt . 



(63) 



The logarithmic divergence of the magnetically broadened bound-bound polarizability 
(see eq. |^0|) removed by the additional Lorentz (or Doppler) broadening. The 
polarizability near the line center can be estimated as 



^,^^(bb) 
AO 



^pafo 



(64) 



for the longitudinal component. A similar estimate is true for the right component, with 
replacing ooq oo+ and /o — >■ /+. 

The above equations completely describe qualitative properties of the cyclic components 
of the polarizability tensor. The results of numerical calculation are shown in Figure 2. 



3.2. Polarization Properties of the Normal Modes 

The NM polarizations can be conveniently described by the parameters q and p defined 
in equations (H) and (12). In wide frequency ranges far from the bound-bound resonances 
and the bound-free thresholds, the hermitian part of the polarizability substantially exceeds 
the antihermitian part, so that |g| 3> \p\, and the polarizations are determined by the 
parameter 

^ = ^~2^' \"^-x"- " • ^^^^ 
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In the case of fully ionized plasma, we have 

At u; 3> uJBp only the electron component is important, and q = Ubc/^ is a smooth function 
of frequency. Kt uj ^ use the parameter g <^ 1, that is, the NMs are circularly polarized in 
a wide range of angles 9. With decreasing frequency, q increases (the polarization ellipses 
become narrower) and reaches its maximum, g^ax — "^^Be/ {'i'^^'^^Bp) — 707 at a; = v^wsp, 
where the NM polarization is virtually linear; then it sharply decreases, crosses zero in the 
very vicinity of proton cyclotron resonance, at ~ ojBpi,^ + '^^Bp/^Be) — 1.0011u;bp and 
tends to — oo (g ~ —ujBe^Bp/^^ at u 0). The fact that q is very small near the proton 
cyclotron resonance means that one should expect a peculiar behavior of the NMs there; 
in particular, strong non-orthogonality of the NMs at some angles. To analyze the NM 
properties around the q = point, one needs to calculate the parameter p whose behavior 
depends on the absorption mechanism operating at these frequencies. 

The relatively simple behavior of q and p for the fully ionized plasma is due to the 
absence of resonances of the longitudinal component of the polarizability xo at finite 
frequencies. For the non-ionized gas, these resonances do exist due to Coulomb binding, 
which strongly complicates the picture (see Fig. 3). Nevertheless, even in this case the NM 
properties can be investigated qualitatively with the aid of approximate expressions for Xa 
derived above. At frequencies well above the photoionization threshold uit ( > 20 — 30 Ry in 
our example), we obtain from equations (pl|), (^), and (|5^) , 



1 = — 

UJ 



1 + 



Cj2 



(67) 



We see that the correction for the longitudinal absorption/refraction increases the value of 
q as compared to the fully ionized case, i. e., the NM polarization is more linear than in the 
fully ionized plasma. At frequencies much lower than ujq — Tm ( ^ 12 Ry in our example), 
and not too close to ti;+, we obtain from equations d^), (|55D, and 



2^ 



fo _^ 3/bf 



(uJo-TMV'io^ 2(30.2-2^2) 



(6^ 



At frequencies much below u^, we have q oc u^/uj, as at high frequencies (eq. |^), but 
this ratio is now multiplied by a small factor ~ fou'^/uQ. Due the Coulomb binding, the 
behavior of q at these low frequencies is qualitatively different from that in the fully ionized 
plasma. The function q{uj) changes its sign near Uj^ and keeps falling with increasing u 
until CO approaches uo — Tm- Equation (B3) is not accurate near u+, where absorption plays 
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a role and the antihermitian component x+i should be taken into account. Nevertheless, it 
satisfactorily describes the qualitative behavior of q even in this region. 

It is clear from equations (|67|) and (^) that q{uj) must change its sign at least once 
more, in a region between ~ (cjq — Ta/) and ~ oot- The behavior of q{uj) in this region 
is determined by interplay of the bound-bound and bound-free refraction and absorption 
(both hermitian and antihermitian components are important), the details depending on 
the broadening mechanisms and the value of the magnetic field. Since q ~ 2xq /{x+i ~ X^i) 
in this region, and x+i ~ X^i (< 0) is a smooth function of frequency, the shape of 
reproduces that of —Xo (cf. Fig. 2). In particular, near the center of an isolated line, 
q = Qmin ^ -0.4/ou;_/rM at u; ~ t^o - ^■4:Tm, and q = g^ax ~ (t^_/2rM) ln(rM/rL) at 
uj ~ ujq. In our example, q reaches its minimum, ~ —237, at hu = 13.35 Ry, then sharply 
grows, changes its sign at hu = 14.22 Ry, reaches the maximum, ~ 472, at hu = 14.42 Ry, 
and then, after some ups and downs associated with the longitudinal bound-free transition, 
approaches the asymptotic curve (|67|). 



The parameter p, which is associated with absorption of radiation (antihermitian 
components), is important whenever \p\ exceeds |g|, that is, in a very narrow range near 

and in the range of strong longitudinal absorption (~ 13 — 20 Ry in our example). 
Around uj+, the antihermitian part is determined by the right-polarization absorption, 
xii > Xo, X-i, and x+i > X-v In this range 

2cosg _2xUxl^X^ 
sin' 9 (x?i)^ + (x^i) 



P — P ^:..2 n — ^ . .H \2 _i_ ( .A ' ^ ^ 



The ratio xii/[{x+i? + xi 

is a smooth function of frequency (the resonances in the 
numerator and denominator cancel each other), and so is Xo ■ In our example, Xo exceeds 
Ix+il in this frequency range (see Fig. 2), and p decreases with decreasing u almost 
smoothly, with a weak wiggle around 5.1 Ry caused by nonmonotonic behavior of the 
relatively small x+i- A characteristic value of p in this range can be roughly estimated as 
p ~ fouj^TM+/{^^o — t^+) (~ 1 in our example — see Fig. 3.). 

In the range of strong longitudinal absorption (and, more generally, at c<j_|_ ^ c<j ^ tu-), 
we can neglect x±i! so that 

p^^l^^-?B}^^, (70) 

X+1 - X 1 ^pa 

and the dependence p(a;) reproduces that of the longitudinal absorption coefficient Hoioj). 
In particular, characteristic values of p near ujq are ~ Jquj^/Tm ~ 5^ (~ 500 in our example). 
They are considerably greater than those near uo^, which results in different behavior of the 
NM polarizations near these frequencies. 



-20- 



In a wide range of frequencies, where |g| ^ |p| ( < 3 Ry, between ~ 6 and ~ 10 Ry, and 
at > 30 Ry, in our example) , the NMs are nearly orthogonal: 

71 T, J- ^ -U 2C0S^ 

" ' 2 ' ' " |g| + (g2 + 1)1/2 = |g| sin2^ + (g2sin^^ + 4cos2^)V2 ' ^ ' 

In the most important domain, a; -C W-, we have |g| ^ 1 in the wide frequency range, 
i. e., the NMs are linearly polarized, iVj] <^ 1, at sin^^ ^ (2/^)^/^. In the ranges of strong 
longitudinal or right-polarization absorption, q and p become comparable, and behavior of 
the NM polarizations is more complicated. The most peculiar behavior occurs near the 
critical frequencies where q = {hujci = 4.751 Ry and hu!c2 = 14.22 Ry in our example). 
For each of these frequencies, there exists an angle 9c such that \p{uJc) I sin^ 9c/ (2 cos 9c) = 1; 
the NM polarizations at a; = a;c and 9 — 9c are completely linear, Vj — 0, and coincide 
with each other. Si — 62 — ±7r/4 (complete non-orthogonahty of NMs). In our example, 
p('^ci) = —3.63 and p{ujc2) = —562, so that the non-orthogonality angles are 9ci — 40.35° 
and 9c2 = 3.41°. Since \p{uOc.,9)\ < 1 dX 9 < 9c., the circular component of polarization 
is present, Vj = ±[(1 — |p|)/(l + bl)]"^'^^, but the orientations of the polarization ellipses 
coincide, 5i = 82 = ±7r/4. If > (i- e., |j9(cjc,^)| > 1)) then the NMs are linearly 
polarized, Vj = 0, but the polarization directions neither coincide nor are orthogonal, 
tan.25j = ±(p2 — l)~^/2. The non-orthogonality points and the curves \q\ = 1 and \p\ = 1 in 
the uj-9 plane are shown in Figure 4. 

In Figure 5 wc demonstrate the frequency dependence of the cllipticitics and position 
angles for three angles of propagation, 9 = 1°, 30°, and 75°, so that wc cover three regions: 
9 < 9c2, 9c2 < 9 < 9ci, and 9 > 9ci- In the case of propagation nearly along the magnetic 
field, the NMs are polarized circularly. The eUipticity Vj decreases shghtly in the range 
around 15 Ry where p is comparable to q. The position angles 5j coincide at the two critical 
frequencies cuci and U!c2- Just above the second critical frequency the position angle varies 
strongly with frequency due to the presence of bound-bound structure in p (see Fig. 3). For 
the intermediate angle 9 = 30°, the polarization of the NMs becomes linear (Vj = 0) at 
the second critical frequency, while at the first one it is elliptical with coinciding positional 
angles. The maximum eUipticity around 5 Ry is due to the fact that in this region p is 
comparable to q, which increases the circular component of polarization. The polarizations 
of the NMs for 9 — 75° > 9ci are close to hnear. They become exactly linear at the two 
critical frequencies cUd and U!c2- Around cUd, however, the polarization of the NMs is 
elliptical, since there is still a significant range where \q\ is small and \p\ is close to unity. 
This is also the cause of the small non-orthogonality of the NMs around ujci- Around U!c2, 
the NMs are polarized linearly and are orthogonal because \p\ is large and the curve q{u!) 
passes the range |g| -C 1 very rapidly. 

The general picture is strikingly different from the case of fully ionized plasma, in 
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which only one critical frequency exists in the vicinity of 00 Bp- The Coulomb binding not 
only shifts this point to higher frequencies and changes the corresponding value of 6c, 
but also creates at least one more critical frequency in the range of strong longitudinal 
absorption. Furthermore, it diminishes \Vj\ at high frequencies, > 30 Ry, and increases it 
at low frequencies, < 5 Ry. 

The location of the critical points where the polarization modes acquire the above- 
described peculiarities is a sensitive function of the temperature and density because it 
depends strongly on the shape of the absorption features. 



3.3. Absorption Coefficients of the Normal Modes 



It follows from the above analysis that in the vicinities of the two points {uJci, ^d) and 
{uJc2, dc2) in the uj-9 plane, the NMs are strongly non-orthogonal (see Fig. 4). Far from these 
non-orthogonality domains, the NM absorption coefficients /ij (j = 1, 2) can be expressed 
as linear combinations of the 'basic' coefficients /Xq, (e. g., Kaminker et al. 1982): 
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As we see from Figure 4, |g| ^ 1 (the NM polarizations are nearly linear) in a wide domain 
of angles and frequencies uj <^ ubc- For this domain we obtain: 



_ fi+i + fi^i fiQ cos^ 9 _ 
2 sm 6 



fi2 = fJ'O sm 6 H cos . 



(75) 



The 'ordinary mode' j = 2 is polarized in the B-q plane. For purely transverse propagation, 
9 = 90°, it is polarized along the magnetic field, so that fj,2 = /io- Since <^ /i+i at 
u <^ UBe (see Fig. 1), the admixture of the transverse-polarization absorption at cos6' 7^ 
is determined by the right-circular absorption coefficient /z+i. The admixture term is 
negligibly small in the domain where equations (^) are applicable, except for a narrow 
frequency range around the resonance of /i+i (~ 2 — 5 Ry in our case). The 'extraordinary 
mode' j = 1 is polarized perpendicular to the magnetic field, and /ii = yU+i/2 at 6* = 90°. 
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Since /io ^ n+i in a wide frequency range, the second term in the equation for /ii, which 
describes the admixture of the longitudinal-polarization absorption, may substantially 
exceed the first term at cos 6' 7^ 0, in spite of the fact that it is proportional to a small factor 



At sin^ <^ q in a wide frequency range, and at sin^ <^ sin^^ near the critical 
frequencies, the NM polarizations are nearly circular (|g| <^ 1 and |p| <^ 1). In this domain. 
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(76) 



(77) 



For purely longitudinal propagation, ^ = 0, the NM polarizations are exactly circular, so 
that fii = and fi2 = fi+i- The left-circular absorption is very weak; in fact, at -C use, 
it is determined by violation of the dipole selection rules due to the motion of atoms in the 
magnetic field (Pavlov & Potekhin 1995), which is neglected in our perturbation approach. 
However, even at very small 6^0, absorption in the j = 1 mode is mainly determined 
by the admixture of the (relatively strong) right-circular and longitudinal absorptions, 
fii ~ (yU_|_i/8) sin"^ 9 + (/io/2) sin^ 9, and the exact value of is not important. 



Equations (|75D, ([76|), and ( |77D allow one to understand qualitatively behavior of 
the spectral and angular dependences of fij depicted in Figures 6, 7, and 8. Figure 6 
demonstrates evolution of fJ^jiuj) with increasing 9 for directions of propagation close to 
the magnetic field. For 9 <^ 9^2 {9 = 1° in Fig. 6), the NMs are polarized almost exactly 
circularly at all the frequencies. The absorption coefficient of the j = 2 mode (dashed curve) 
is determined by the right-circular absorption /i+i (dominates a.t hu <8 Ry [bound-bound 
transition to (010)] and > 15 Ry [bound-free transition to (01e-|-)] — cf. Fig. 1) and by a 
small addition, (/io/2) sin^ 9 of the longitudinal absorption (dominates at the longitudinal 
resonance, hu ^ 10 — 15 Ry). The other mode is solely determined (at this angle) by the 
same contribution of the longitudinal absorption. The absorption coefficients fii and fi2 
almost coincide with each other at the region of strong longitudinal absorption. 

The increase of 9 towards and beyond the critical value 9c2 = 3.41° is accompanied 
by fast growth of both NM absorption coefficients in the region of prevailing longitudinal 
absorption, hu > 10 Ry; the peak at /icj ~ 14 — 15 Ry becomes ~ 40 times higher when 9 
changes from 1° to 5°. The coefficients also approach each other at these frequencies when 
9 approaches 9^2, and separate again when 9 moves away from 9c2 (see also Fig. 8). The 
resonance of right-circular absorption around 5 Ry remains almost unchanged in one mode, 
and its admixture to the other mode grows oc sin^ 9. 
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Figure 7 demonstrates behavior of Hj{uj) at intermediate 6. The left panel {6 = 30°) 
shows a typical frequency dependence at 6c2 < < dci- The NM polarization is almost 
circular around the critical frequency hud = 4.75 Ry, and the absorption coefficients can 
be evaluated there from equations (7B) and (|77p. The main distinction from the case of 



small 6 is that the right-circular absorption gives a larger contribution to /ii, so that the 
resonance in that NM mode became much more pronounced. At hu >6 Ry and < 4 Ry, 
the NMs are polarized almost linearly, (see Figure 5), and their absorption coefficients are 
given by equation (^). It should be noted, however, that fi2 in equation ([75|) corresponds 
to the solid line aX hu >6 Ry, and to the dashed line at hu <A Ry. The reason is that in 
the presence of the critical frequencies and non-orthogonality points a unique labeling of 
the NMs along the continuous curves of the NM absorption and refraction is impossible (see 
Pavlov et al. 1980); this fact, however, does not impose any difficulties because both NMs 
participate in the radiative transfer, and their labeling is purely conventional. 

The angle 6 = 41° only slightly exceeds 6ci = 40.35°. At this angle the NMs are 
substantially non-orthogonal near Ud, and the absorption coefficients nearly coincide in this 
region of the right-circular resonance (the middle panel of Fig. 7). At frequencies below and 
above this region, the NM polarization is close to linear, and the NM behavior is similar to 
that for ^ = 30°. 

At the angle 6 = 75°, which is well above 6ci, the NM polarizations are almost linear at 
all frequencies, although there is a very narrow region of enhanced ellipticity at /icu ^ 5.1 Ry 
(see Fig. 5). The absorption coefficients are fairly accurately described by equation ( [75D {fi2 
corresponds to the solid curve for all the frequencies). Notice that /ii and fi2 coincide just 
in the enhanced ellipticity region. 

Figure 8 presents the absorption coefficients function of the propagation angle 

for a few characteristic frequencies. At low frequencies u < Ud (see the panel for hu = 2 Ry 
in Fig. 8), where the absorption coefficients are determined by the magnetically broadened 
wings of the bound-bound lines, the NM polarization changes from circular to linear with 
increasing angle. Absorption in one mode slightly decreases with 6, while in the other mode 
it grows very fast at small angles because of increasing contribution of /iq (see eq. [ffTf); 



the coefficients become close to each other at ^ > 20°, where the NMs are almost linearly 
polarized. 

The picture changes considerably when we approach the first critical frequency, 
huci = 4.751 Ry. Here, for 6 = Oci = 40.35°, the (linear) polarizations and the absorption 
coefficients of the modes coincide with each other: 6i = 62 = tt/A and Vi = V2 = 0. The 
modes are linearly polarized at angles above ^ci, and they become more orthogonal with 
increasing 9. Below Od-, their polarizations are elliptical, and the ellipticity increases with 
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decreasing angle. 

In the region between the critical frequencies, Ud < uo < uj^i^ = 7 Ry in Fig. 8), 
the longitudinal absorption fiQ (the low-frequency wing of the longitudinal resonance) 
is much greater than the right-circular absorption /x+i (the high-frequency wing of the 
spectral line and the low-frequency wing of the photoionization edge). The NM polarization 
is nearly circular at small 6 ( < 10° for huj = 7 Ry), and, in accordance with equations 
([f6|) and (ffTl), both /ii and fj,2 grow with 9: ni ~ /i+i + (/Uo/2) sin^ ^, /i2 — (/io/2) sin^ 6*. 
The polarization becomes almost linear at > 20°, where /i2 — /iosin^^ keeps growing, 
whereas /ii = (/i+i/2) + (/ig/g^) cot^ 6* decreases because of diminishing contribution from 
the longitudinal absorption. 

At the second critical frequency TiuJc2 = 14.22 Ry, the NMs coincide at the critical angle 
9c2 = 3.41°. For < 6c2, the NMs modes are nearly circular, and the absorption coefficients 
increase with angle, being very close to each other (cf. Fig. 6) because the both are mainly 
contributed from the longitudinal absorption (/xo/z^+i — 1-3 x 10^ at this frequency). For 
6 > 6c2, the linearly polarized NMs become more orthogonal with increasing 6: ^2 ^ and 
6i 7r/2; absorption coefficients of the ordinary and extraordinary modes grow up to /iq 
and decrease down to /i+i/2, respectively. 

Behavior of fij{6) above the second critical frequency is demonstrated in Figure 8 
for hu = 14.82 Ry, near the peak of longitudinal absorption, and for hu = 30 Ry, where 
both iiq and /xi are determined by the bound-free transitions. At these frequencies, the 
NM polarizations are virtually orthogonal for any 6. Since |g| ^ 1, the polarizations 
are circular in a narrow range of angles, sin^ <^^^^^, and they are linear in a wide 
range of angles. The absorption coefficient of the ordinary mode grows monotonically 
with 9: fi2 — fJ^+i + (/io/2) sin^ 6' at small 9, and /i2 — /iosin^6' at larger angles. For the 
extraordinary mode, /ii ~ (/Uo/2) sin^ 9 grows at small 9, reaches its maximum value, ~ /io/*? 
at sin 6* ~ and then decreases gradually, /ii ~ (/io/*?^) cot^ 9 + /i+i/2. 



4. DISCUSSION 



In this paper we suggested a convenient method for calculating the absorption 
coefficients of the normal modes in a magnetized medium. Given the 'basic' absorption 
coefficients fiQ (for the linear polarization along B at transverse propagation) and (for 
the circular polarizations at longitudinal propagation), we obtain the antihermitian part 
of the polarizability tensor and use the Kramers-Kronig transformation to calculate its 
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hermitian part. The knowledge of the full polarizability tensor allows us to calculate the 
polarizations and the absorption coefficients of the NMs for any frequency and direction of 
propagation. We applied the method to a nonionized, strongly magnetized hydrogen gas 
for conditions expected in atmospheres of neutron stars. Absorption and propagation of 
radiation in such a medium is determined by the bound-bound and bound-free transitions. 
The corresponding spectral lines and photoionization edges are strongly broadened due to 
distortion of the structure of atoms moving across the magnetic fields, and the shapes of 
the lines depend essentially on the direction of propagation and polarization of radiation. 
As a result, the angular and spectral dependences of the NM characteristics are much 
more complicated than in a fully ionized plasma. In particular, the presence of the 
nonionized component leads to appearance of additional critical frequencies and angles, 
where peculiarities of the NM polarizations cause additional crossings of the curves /ijico) 
or fij{9). 

The resulting absorption coefficients can be used to solve the radiative transfer 
equations for the NMs in neutron star atmospheres and to calculate spectra, angular 
distribution and polarization of emergent radiation. Since the temperature grows inward in 
the atmosphere, the spectrum and direction of the radiation emitted is mainly determined 
by the NM which has lower opacity and, therefore, escapes from deeper and hotter layers. 
We see from Figure 6 that at small 9 the NM of lower opacity shows a strong, asymmetric 
resonance associated with the bound-bound and bound-free transitions for the longitudinal 
polarization. Consequently, we should expect a broad, asymmetric absorption feature (at 
hcu 14 — 15 Ry in our example) in the emergent spectrum observed at directions close 
to the magnetic field direction (c. g., when the neutron star is seen pole-on). Figure 7 
shows that at intermediate and large 9 the strongest absorption in the NM of lower opacity 
is that at the bound-bound resonance for the right-circular polarization (at ~ 5 Ry), so 
that the corresponding absorption feature should be most clearly seen at these directions. 
From the angular dependences of nj shown in Figure 8 we see that at most frequencies 
the NM of lower opacity sharply grows with 9 at small 9, reaches a maximum, and then 
decreases gradually when 9 — > 90°. This means that the angular distribution I{9) of the 
intensity of radiation is beamed along the magnetic field (a 'pencil component' of the 
angular distribution); a typical width of the beam is < g"^/^ at g » 1 (A^ ~ 3° — 10° at 
hu ~ 6 — 30 Ry in our example). The maximum of /Xj(^) should correspond to a minimum 
of I{9), i. e., a broad 'fan component' in the plane perpendicular to the magnetic field is 
emitted in addition to the sharp 'pencil component'. The pencil component broadens when 
the frequency approaches ujcu at a; < ujci the minimum of I{9) (and the fan component) 
disappears. The polarization of the emergent radiation is expected to be strongly polarized 
at frequencies and angles where the the NM absorption coefficients are substantially 
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different. The polarization is circular in the narrow pencil component, and it is linear 
(perpendicular to the magnetic field) in the fan component. 

In the present paper we considered only the case of low temperatures, when the atoms 
are expected to be nonionized and the effect of their motion on the radiative transitions 
can be considered in the framework of the perturbation theory. At higher temperatures, 
the ionized component of the plasma would contribute to the polarization tensor, NM 
polarizations, and NM opacities. This will result in even more complicated spectral and 
angular behavior of the NM characteristics. To include both ionized and nonionized 
components, one should first solve the nontrivial problem of ionization equilibrium of 
a strongly magnetized, highly nonideal plasma. Including the (nonperturbative) effects 
of the atomic motion at higher temperatures (Bezchastnov & Potekhin 1994; Pavlov & 
Potekhin 1995) would be also desirable, particularly for propagation at directions close to 
the magnetic field. Besides, it is possible that not only hydrogen, but also other atoms 
(e.g. helium, or even iron) can be present in neutron star atmospheres. The study of the 
structure and radiative transitions in the heavier atoms in very strong magnetic fields is a 
very complicated problem. When these problems are solved, the results can be easily used 
to calculate the NM properties with the method suggested in the present paper. 

The research has been partially supported through NASA grants NAG5-2807, 
NAG5-2868 and NASW-4690. We are grateful to Alexander Potekhin who supplied the 
codes for calculating the photoionization cross section used in this paper. We thank Victor 
Bezchastnov, Cole Miller, and Yura Shibanov for useful discussions. 
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Table 1: Bound-bound and bound- free transitions from the ground state, their oscilator 
strengths, and magnetic, Lorcntz, and Doppler widths of the associated lines and photo- 
ionization edges at /cT = 1 Ry, p = 1 g cm~^, and B = 2.35 x 10^^ G. 
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Figure Captions 

Figure 1. Frequency dependences of the absorption coefficients /Xq, for the three 
polarizations: a — (hnear polarization along the magnetic field B), and a — ±1 
(circular polarizations in the plane perpendicular to B). The right panel shows detailed 
structure of /Iq and in the range of bound-bound resonances and bound-free edges. 
The contributions of the bound-bound (bb) and bound-free (bf) transitions are shown 
separately. The low-frequency wings of the resonances are determined by the magnetically 
broadened spectral lines. 

Figure 2. The upper panel shows the circular components x±i of the hermitian and 
antihermitian parts polarizability tensor. The left-circular component is very small in 
this frequency range. The major features of x+ and x+ are associated with the bound-bound 
resonance at 5.1 Ry. The lower panel shows Xo and Xo ■ "^^^ main features (shown in 
the inset separately) are due to the photoionization edge at 15.3 Ry and bound-bound 
resonance at 14.4 Ry. 

Figure 3. The upper panel shows the global behavior of the parameters q and p which 
determine the NM polarization. The two lower panels show q and p in the vicinity of the 
critical frequencies where q = 0. Note the different vertical scales. 

Figure 4. The solid and dashed curves connect the points in the uj-6 plane for which 
|g| = 1 and \p\ = 1, respectively. The labels ci and C2 denote the non-orthogonality points 
where q = and \p\ = 1. 

Figure 5. The left three panels show ellipticity and the right three panels show position 
angles of the NM polarization ellipses as functions of frequency for three angles between B 
and direction of propagation: 9 = 1°, 30°, and 75°. Note the different scale in the lower 
panels. 

Figure 6. Absorption coefficients of the NMs vs. frequency for three directions of 
propagation close to B. 

Figure 7. Absorption coefficients of the NMs vs. frequency for three intermediate 
angles 9. 

Figure 8. Angular dependences of the absorption coefficients of the NMs for six photon 
energies. 
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